**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



Spectroscopic Analysis of subluminous B Stars in Binaries 
with compact Companions 



S. Geier 1 , C. Karl 1 , H. Edelmann 2 , U. Heber 1 and R. Napiwotzki 3 

. 1 Dr.-Remeis-Sternwarte, Institute for Astronomy, University 

Erlang en- Nuremberg, Sternwartstr. 7, 96049 Bamberg, Germany 
2 McDonald Observatory, University of Texas at Austin, 1 University 
Station, C1402, Austin, TX 78712-0259, USA 



Q . 3 Centre of Astrophysics Research, University of Hertfordshire, College 

q ; Lane, Hatfield ALIO 9 AB, UK 

Abstract. The masses of compact objects like white dwarfs, neutron stars 
and black holes are fundamental to astrophysics, but very difficult to measure. 
We present the results of an analysis of subluminous B (sdB) stars in close 
binary systems with unseen compact companions to derive their masses and 
■ clarify their nature. Radial velocity curves were obtained from time resolved 

^C) \ spectroscopy. The atmospheric parameters were determined in a quantitative 

. spectral analysis. Based on high resolution spectra we were able to measure the 

fN) ' projected rotational velocity of the stars with high accuracy. The assumption 

of orbital synchronization makes it possible to constrain inclination angle and 
companion mass of the binaries. Five invisible companions have masses that 
are compatible with that of normal white dwarfs or late type main sequence 
^ (— i stars. But four sdBs have very massive companions like heavy white dwarfs 

O (> 1M©), neutron stars or even black holes. Such a high fraction of massive 

I compact companions is not expected from current models of binary evolution. 

CO 

1. Introduction 

The mass of a star is it's most fundamental parameter. However, a direct mea- 
surement is possible in some binary stars only. Eclipsing, double lined systems 
are first choice. White dwarfs, neutron stars and stellar black holes are the 
aftermath of stellar evolution. In binaries such faint, compact objects are out- 
shined by their bright companions and therefore their orbital motion cannot be 
measured. As a consequence only lower limits to the companion mass can be 
derived. With the analysis method shown here, these limitations can partly be 
overcome. 

Subluminous B stars, which are also known as hot sudwarf stars (sdBs), show the 
same spectral characteristics as main sequence stars of spectral type B, but are 
much less luminous. They are considered to be helium core burning stars with 
very thin hydrogen envelopes and masses around 0.5 Mq. Different formation 
channels have been discussed. As it tur ned out, a large frac tion of the sdB stars 
are members of short period binaries ( Maxted et al.l 20011 ). For the se systems 



comm on envelope ejection is the most probable formation channel (jHan et al 
2002). In this scenario two main sequence stars of different masses evolve in 
a binary system. The heavier one will first reach the red giant phase and fill 
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its Roche lobe. If the mass transfer to the companion is dynamically unsta- 
ble, a common envelope is formed. Due to friction the two stellar cores loose 
orbital energy, which is deposited within the envelope and leads to a shortage 
of the binary period. Eventually the common envelope is ejected and a close 
binary system is formed, which contains a helium core burning sdB and a main 
sequence companion. When the latter reaches the red giant branch, another 
common envelope phase is possible and can lead to a close binary with a white 
dwarf companion and an sdB. All known companions of sdBs in such systems 
are white dwarfs or late type main sequence stars. If massive stars are involved, 
the primary may evolve into a neutron star (NS) or a black hole (BH) rather 
than a white dwarf. Since massive stars are very rare, only few sdB+NS or 
sdB+BH systems are expected to be found. 

Since the spectra of the program stars are single-lined, they reveal no informa- 
tion about the orbital motion of the sdBs' companions, and thus only their mass 
functions can be calculated. 

M c 3 omp sin 3 i = PK 3 
/m (M comp + M sdB ) 2 2nG • 1 > 

Although the RV semi-amplitude K and the period P are determined by 
the RV curve, M S( jb ; Mm rnp and s'mi rem ain free parameters. Binary pop- 
ulation synthesis models (jHan et al.1 120021 ) indicate a possible mass range of 
-^sdB = 0.30—0.48 M for sdBs in binaries, which underwent the common en- 
velope ejection channel. The mass distribution shows a sharp peak at about 
0.46 M Q . 

For close binary systems, the components' stellar rotational velocities are con- 
sidered to be tidally locked to their orbital motions, which means that the or- 
bital period of the system equals the rotational period of the companions. If 
the companions are synchronized in this way the rotational velocity i> rot can be 
calculated. 

v m t = — p — • (2) 
The stellar radius R is given by the mass radius relation. 



R = sj^f- (3) 

The measurement of the projected rotational velocities v ro t sin i therefore 
allows to constrain the systems' inclination angles i. With M sc jb as free param- 
eter the mass function can be solved and the inclination angle as well as the 
companion mass can be derived. Because of sini < 1 a lower limit for the sdB 
mass is given. To constrain the system parameters in this way it is necessary to 
measure K, P, log <? and u rot sini with high accuracy. 



2. Observations and Radial Velocity Curves 



Ten stars were observed at least twice with the high resolution spectrograph 
UVES at the ESO VLT. Additional observations were made at the ESO NTT 
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(equipped with EMMI), the Calar Alto Observatory 3.5 m telescope (TWIN) and 
the 4 m WHT (ISIS) at La Palma. Two of the stars (PG 1232-136, TONS 183) 
were observed with the high resolution FEROS instrument at the 2.2 m ESO 
telescope at La Silla. 

Radial velocities of the individual observations were determined by calculat- 
ing the shifts of the measured wavelengths relative to their laboratory values. 
Results are displayed in Table 1. 



3. Atmospheric parameters and projected rotational velocities 

The spectra were corrected for the measured RV and coadded. Atmospheric 
parameters were determined by fitting simultaneously each observed hydrogen 
and helium line with a grid of metal-line blanketed LTE model spectra. Results 
are listed in Table 1. 



Table 1. Effective temperatures, surface gravities, orbital periods, radial ve- 
locity semi-amplitudes and projected rotational velocities of the visible com- 
ponents. Typical error margins for T e ff and logg are 500 K and 0.05 dex 
respectively, f Orbital parameters of this system are still preliminary. 



System 






logff 


P 


K 


v m t sin i 






[K] 




[d] 


[km s _1 ] 


[kms -1 ] 


HE 0532- 


-4503 


25390 


5.32 


0.26560 ± 0.00010 


101.5 ± 0.2 


11.1 ± 0.6 


PG1232- 


-136 


27500 


5.62 


0.36300 ± 0.00030 


129.6 ± 0.04 


6.2 ± 0.8 


WD 0107 


-342t 


24300 


5.32 


0.38000 ± 0.00050 


135.0 ± 1.0 


20.4 ± 0.9 


HE 0929- 


-0424 


29470 


5.71 


0.44000 ± 0.00020 


114.3 ± 1.4 


7.1 ± 1.0 


HE0230- 


-4323 


31100 


5.60 


0.44300 ± 0.00050 


64.1 ± 1.5 


12.7 ± 0.7 


TONS 183 


26100 


5.20 


0.82770 ± 0.00020 


84.8 ± 1.0 


6.7 ± 0.7 


HE 2135- 


-3749 


30000 


5.84 


0.92400 ± 0.00030 


90.5 ± 0.6 


6.9 ± 0.5 


HE 1421- 


-1206 


29570 


5.55 


1.18800 ± 0.00100 


55.5 ± 2.0 


6.7 ± 1.1 


HE 1047- 


-0436 


30242 


5.66 


1.21325 ± 0.00001 


94.0 ± 3.0 


6.2 ± 0.6 


HE2150- 


-0238 


30200 


5.83 


1.32090 ± 0.00500 


96.3 ± 1.4 


8.3 ± 1.3 


HE 1448- 


-0510 


34690 


5.59 


7.15880 ± 0.01300 


53.7 ± 1.1 


6.7 ± 2.5 


WD 0048 


-202 


29960 


5.50 


7.44360 ± 0.01500 


47.9 ± 0.4 


7.2 ± 1.3 



In order to derive v TO t sin i, we compared the observed spectra with rota- 
tionally broadened, synthetic line profiles. The latter ones were computed using 
the LINFOR program. Since sharp metal lines are much more sensitive to ro- 
tational broadening than Balmer or helium lines, all visible metal lines were 
included. A simultaneous fit of elemental abundance and projected rotational 
velocity was performed separat ely for every identified line using the FITSB2 
routine ( Napiwotzki et al. 20041 ). Mean value and standard deviation were cal- 



culated from all measurements. Seeing induced variations in the instrumental 
profile and the noise level were the dominant error sources. Information on the 
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Figure 1. Companion mass as a function of primary (sdB) mass of the bi- 
naries HE 0532-4503 (left hand panel) and PG 1232-136 (right hand panel). 
The horizontal line marks the Chandrasekhar limit. The dotted vertical lines 
mark the theoretica l sdB mass range for the common envelope ejection chan- 
nel (jHan et alJl2002f ). 



actual seeing conditions for every exposure have been extracted from the ESO 
seeing monitor archive. All other possible sources of systematic errors turned 
out to be negligible. 



4. Nature of the unseen Companions 

Knowing P, K, log g and v ro t sin i we can calculate the mass of the unseen com- 
panion from equations (1) - (3) for any giv en primary mass (see Fig. [1]). We 
adopt the mass range from ( Han et al. 20021 ). marked by the dotted lines in Fig. 



[TJ To constrain the comp anion masses and inclination angles we adopted the 
mass range for sdBs from Han et al. ( 20021 ). There are no spectral signatures 



of companions visible. Main sequence stars with masses higher than 0.45 Mq 
could therefore be excluded because of their high luminosities in comparison to 
the sdB stars. The possible companion masses can be seen in Tab. 2. Four of 
the analysed systems have companion masses, which are compatible with either 
typical white dwarfs (WD) or late main sequence stars (late MS). The compan- 
ion of WD 0107—0342 is a white dwarf of intermediate mass. 
The very similar HE 0929—0424 and TONS 183 have to have quite massive 
companions. If their sdB primaries are of low mass (~ 0.3 Mq) the compan- 
ions could be a heavy white dwarfs. At the most probable sdB mass, how- 
ever, their companions would exceed the Chandrasekhar mass limit. There are 
only two kinds of objects known with such high masses and such low luminosi- 
ties - neutron stars (NS) and stellar mass black holes (BH). The two systems 
HE 0532— 4503 and PG 1232— 136 have even higher companion masses, which 
would exceed the Chandrasekhar limit in any case (see Fig. [I]) . The three bina- 
ries HE 1448-0510, HE 2150-0238 and WD 0048-202 could not be solved with 
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the described method. The minimum sdB masses exceeded 1.3 Mq and were not 
consistent with the theoretical mass range. These systems cannot be synchro- 
nized. 

Binaries hosting a neutron star or a black hole are a very rare class of objects. 
From about 50 analysed sdB binaries in our samples, we found two, possibly four 
candidate systems. This fraction of 4-8 % is much too high to be compatible 
with any binary evolution model known so far (Podsiadlowski priv. comm.). 
The reliability of our method has to be discussed again. Only the high preci- 
sion of the v m t sin i measurements made it possible to constrain the parameters 
reasonably. The projected rotational velocities were low and very close to the 
detection limit. As described above we tried to quantify all possible system- 
atic effects and the overall results were very consistent. Only slight changes in 
v r otsini would lead to inconsistent solutions, because the method is very sen- 
sitive to this parameter. But even if there would be unaccounted systematic 
effects (e.g. short period pulsations), they would always cause an extra broad- 
ening of the lines. The measured broadening is then due to rotation and the 
unaccounted effects, which means the deduced rotational broadening is higher 
than the real one. All systematic effects lead to lower v rot sini. But the un- 
expectedly high masses of the companions are caused by the unexpectedly low 
measured projected rotational velocities as can be seen from the equations (1) 
- (3). Unaccounted systematic effects would therefore lead to even higher com- 
panion masses. 

Our method rests on the assumption of orbital synchronization. Many eclipsing 
binaries in close orbits were observed to be synchronized. At least two short 
period (~ 0.1 d) sdB binaries with white dwarf companions are proven to be 
synchronized ( Orosz Sz Wadd[l999l . Geier et al. these proceedings). 
The question, which mechanism is responsible for this effect, is not yet set - 
tled . The approximation form ulas for the synchronization time of IZahn ( 19771 ) 
and lTassoul Tassoull (1992) for stars with radiative envelopes and convective 
cores were used to compare with the lifetime on the Extreme Horizontal Branch 
(EHB) £ehb ~ 10 8 yr. Nine short period binaries can be solved consistently and 
have synchronization times much lower than their EHB lifetimes. Three systems 
could not be solved consistently. These three binaries have the longest orbital 
periods and synchronization times near or above the EHB lifetime. The two 
systems with the most massive companions have the shortest orbital periods 
and therefore should most likely be synchronized (see Table 1 & 2). 
But as long as the question of tidal synchronization is not settled (see review by 



Zahnll2005l ). all timescales have to be taken with caution. Detailed calculations, 



which take the internal structure of sdBs into account, are not available and ur- 
gently needed. From our results we draw the conclusion, that the assumption of 
tidally locked rotation is reasonable in the case of close sdB binaries for orbital 
periods shorter than 1.3 d. 



5. Conclusion 

Out of 12 analysed sdB binaries, five have companion masses compatible with 
white dwarfs of typical masses or late type main sequence stars. These systems 
are in full agreement with binary population synthesis simulations. Four bina- 
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ries have surprisingly high companion masses, which leads to the conclusion, 
that the companions have to be white dwarfs of unusually high mass or even 
neutron stars or black holes. This high fraction cannot be explained with cur- 
rent evolutionary calculation. As our analysis assumes synchronization a better 
understanding of this process for sdB stars is urgently needed. More observa- 
tions of close sdB binaries and candidate systems with high and low resolution 
spectrographs are needed to enhance the size of our sample. 
The presence of such a high fraction of heavy binaries in our samples raises 
several questions. Is the formation and evolution of sdB stars linked to that of 
heavy compact objects like neutron stars or black holes? Can sdB stars be used 
as tracers to find more of these exotic objects? Is there a hidden population of 
these objects, which is not visible through strong X-ray emission? 



Table 2. Inclination angles, rotational velocities, companion masses and 
possible nature of the unseen companions. The lower companion mass corre- 
sponds to an sdB of 0.3 Mq, the higher limit to an sdB of 0.48 M f Derived 
parameters of this system are still preliminary. 



System 




i 

[de { 




Vrot 

[kms -1 ] 


■'"comp 
[Mq] 


Companion 


HE 0532- 


-4503 


13 - 


17 


47 


1.40 - 


3.60 


NS/BH 


PG1232- 


-136 


14 - 


19 


25 


2.00 - 


7.00 


NS/BH 


WD 0107 


-342t 


17 - 


22 


33 


0.45 - 


0.95 


WD 


HE 0929- 


-0424 


23 - 


29 


18 


0.60 - 


2.40 


WD/NS/BH 


HE0230- 


-4323 


38 - 


50 


21 


0.18 - 


0.35 


WD/late MS 


TONS 183 


22 - 


29 


18 


0.60 - 


2.40 


WD/NS/BH 


HE 2135- 


-3749 


66 - 


90 


8 


0.35 - 


0.45 


WD/late MS 


HE 1421- 


-1206 


56 - 


90 


8 


0.15 - 


0.30 


WD/late MS 


HE 1047- 


-0436 


62 - 


90 


7 


0.35 - 


0.60 


WD/late MS 


HE2150- 


-0238 












no solution 


HE 1448- 


-0510 












no solution 


WD 0048 


-202 












no solution 
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